function [X_norm, mu, sigma] = featureNormalize(X)
    m = size(X, 1);
    mu = mean(X);
    sigma = std(X);
    X_norm = (X - repmat(mu, m, 1)) ./ repmat(sigma, m, 1);
end 